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This  report  consists  of  two  parts,  tn  Part  I*  a  time  dependent  form  of 
the  reduced  wave  equation  of  Berkhoff  is  developed  for  the  case  of,wa*ve« 
propagating  over  a  bed  consisting  of  ripples  superimposed  on  an  otherwise 
slowly  varying  mean  depth  which  satisfies  the  mild  slope  assumption.  The 
rlppl*s  are  assumed  to  have  wavelengths  on  the  order  of  the  surface  wave 
length  but  amplitudes  which  scale  as  a  small  parameter  along  with  the  bottom 
slope.  The  theory  is  verified  by  showing  that  It  reduces  to  the  case  of  plane 
waves  propagating  over  a  one -dimensional ,  infinite  patch  of  sinusoidal 
ripples,  studied  recently  by  Davies  and  Heathershaw  and  Met.  we  then  study 
two  cases  of  interest*^  formulation  and  use  of  the  coupled  parabolic  equations 
for  propagation  over  patches  of  arbitrary  form  In  ord^r  to  study  wave 
reflection,  and  propagation  of  trapped  waves  along  an  infinite  ripple  patch. 

In  the  second  part,  we  use  the  results  of  Part  1  to  estend  the  results 
for  weakly-nonlinear  wave  propagation  to  the  case  of  partial  reflection  from 
bottoms  with  mild-sloping  mean  depth  with  superposed  small  amplitude 
undulations. 
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ABSTRACT 


This  report  consists  of  two  parts.  In  Pert  l,  e  t tee  dependent  fore  of 
the  reduced  wave  equation  of  Berkhoff  Is  developed  for  the  case  of  waves 
propagating  over  a  bed  consisting  of  ripples  superimposed  on  an  otherwise 
slowly  varying  mean  depth  which  satisfies  the  at  Id  elope  assumption.  The 
ripples  are  assumed  to  have  wavelengths  on  the  order  of  the  surface  wave 
length  but  amplitudes  which  scale  as  a  small  parameter  along  with  the  bottom 
slope.  The  theory  Is  verified  by  showing  that  It  reduces  to  the  case  of  plane 
waves  propagating  over  a  one-dlmenslonal,  infinite  patch  of  sinusoidal 
ripples,  studied  recently  by  Davies  and  Heatherahaw  and  Met.  Me  then  study 
two  cases  of  Interest;  formutation  and  use  of  the  coupled  parabolic  equations 
for  propagation  over  patches  of  arbitrary  fora  In  order  to  study  wave 
reflection,  and  propagation  of  trapped  waves  along  an  Infinite  ripple  patch. 

In  the  second  part,  we  use  the  results  of  Part  I  to  extend  the  results 
for  weakly -non linear  wave  propagation  to  the  caae  of  partial  reflection  fros 
bottoms  with  mlld-sloplng  mean  depth  with  superposed  small  amplitude 
undulations. 
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Chapter  I 


Sumraa ry 

Wave  models  for  waves  propagating  over  uneven  topography  have  undergone  a 
great  deal  of  development  In  the  recent  past,  with  the  ma^rity  of  work 
centered  on  the  extreme  examples  of  very  abrupt  topography,  where  the  bottom 
may  be  reasonably  schematized  by  sections  of  constant  depth  separated  by 
discontinuities,  and  of  topography  with  very  mild  bottom  slopes,  where  the 
depth  Is  assumed  to  change  very  slowly  over  the  space  of  a  wavelength. 

Although  some  of  the  work  encompassed  by  the  ongoing  studies  under  the  present 
ONR  support  have  dealt  with  the  first  example  (Kirby  and  Dalrymple,  1983a, 
plus  continuing  work  on  trenches  and  dredged  channels  with  currents),  most  of 
the  effort  In  the  present  contract  period  has  dealt  with  the  mlld-slope 
extreme,  where  the  basic  formulations  are  either  In  the  form  of  reduced  wave 
equations  (Berkhoff,  1972;  Kirby,  1984)  or  WKB  formulations  giving  evolution 
equations  for  the  slowly  varying  domain  (Chu  and  Mei,  1970;  Djordjevlc  and 
Redekopp,  1978). 

This  report  presents  two  theoretical  results  which  extend  the 
capabilities  for  modelling  waves  in  shallow  water.  Chapter  II  presents  a 
linear,  reduced  wave  equation  which  Is  applicable  to  the  case  of  waves 
propagating  over  bottom  undulations  which  are  small  in  amplitude  but  which  may 
have  length  scales  on  the  order  of  the  surface  wavelength.  This  Eormulation 
extends  the  basic  mild-slope  equation  to  the  case  of  abrupt  bars  or  similar 
topography  resting  on  an  otherwise  mildly  sloping  bottom.  Comparisons  to 
previous  theoretical  and  experimental  results  are  given. 

In  Chapter  III,  the  results  of  Chapter  II  are  combined  together  with  the 
coupled  forward-  and  backward-scattered  wave  formulation  of  Liu  and  Tsay 


(1983)  and  the  perturbation  method  based  on  a  Lagrangian  formulation  (Kirby, 
1983)  in  order  to  develop  a  model  for  the  gradual  reflection  of  weakly 
nonlinear  Stokes  waves.  Numerical  results  are  given  for  several  examples 
involving  waves  propagating  in  one  and  two  dimensions. 
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Chapter  II 


A  General  Wave  Equation  for  Waves  over  Rippled  Beds 

. I  Introduction 

The  problem  of  reflection  of  surface  waves  by  patches  of  large  bottom 
dulations  has  received  an  Increasing  amount  of  attention  recently,  due  to 
is  mechanism's  possible  Importance  In  the  development  of  shore-parallel 
rs.  Davies  and  Heathershaw  (1984)  and  Mel  (1984)  have  recently  studied  the 
ise  of  reflection  from  sinusoidal  topography  and  have  provided  analytic 
eatments  which  elucidate  the  mechanism  of  a  resonant  Bragg  reflection  at  the 
>int  where  the  wavelength  of  the  bottom  undulation  is  one  half  the  wavelength 
:  the  surface  wave.  The  analytic  results  of  both  studies  are  seen  to  agree 
iite  well  with  the  laboratory  data  of  Davies  and  Heathershaw  (1984, 
riginally  presented  by  Heathershaw,  1982).  Based  on  the  analytic  and 
icperimental  results,  Mei  suggests  that  reflection  of  waves  from  an  initially 
solated  bar  (such  as  a  break-point  bar  formed  at  the  outer  edge  of  the 
arfzone)  may  be  sufficient  to  induce  the  formation  of  bars  further  offshore 
f  the  initial  bar. 

The  analytic  results  presented  to  date  illustrate  the  major  features  to 
e  expected  when  studying  reflection  from  a  system  of  bars.  However,  they  are 
oo  limited  in  scope  to  provide  a  direct  treatment  in  the  case  of  natural  bed 
orms  varying  arbitrarily  in  two  horizontal  directions.  For  this  reason,  the 
resent  study  concentrates  on  the  development  of  a  general  wave  equation  which 
s  applicable  to  Linear  surface  waves  in  Intermediate  or  shallow  water 
epths.  The  resulting  equation  is  similar  in  spirit  to  the  reduced  wave 
quation  of  Berkhnff  (1972),  but  extends  the  usual  mild-slope  approximation  to 
nclude  rapidly-varying,  small-amplitude  deviations  from  the  slowly-varying 


ean  depth. 
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After  deriving  the  general  equation  in  Section  2,  we  consider  the 
;orrespondence  between  the  present  results  and  those  of  Mel  (1984)  in 
Section  3.  A  simple  numerical  scheme  is  then  used  to  solve  the  reduced  wave 
aquation  for  the  case  of  reflection  over  one-dimensional  topographies. 

In  Section  4,  we  apply  a  splitting  method  in  order  to  reduce  the  elliptic 
form  to  two  coupled  parabolic  equations  for  forward  and  backseat i ered  waves. 
The  resulting  equations  for  amplitude  of  the  forward  and  backward  propagating 
waves  extend  the  results  of  Mel  to  the  case  of  arbitrary  topographic 
variations  and  include  possible  diffraction  effects. 

Finally,  in  Section  5  we  consider  the  possibly  Interesting  rase  of 
trapped  waves  propagating  along  a  ripple  patch.  The  problem  Is  developed  is 
an  eigenvalue  problem  for  the  wavelength  along  the  patch,  which  Is  then  solved 
numerically  for  the  case  of  arbitrary  topography,  following  the  method  of 
Kirby,  Dalrymple  and  Liu  (1981). 

II. 2  Derivation  of  the  Wave  Equation 

The  depth-integrated  wave  equation  for  monochromatic,  linear  waves 
propagating  over  small-amplitude  bed  undulations  may  be  formulated  following 
either  a  variant  of  the  Green's  formula  method  of  Smith  and  Sprlnks  (1975)  or 
Liu  (1983),  or  by  using  the  Lagranglan  formulation  of  Kirby  (1984).  Here  the 
Green's  formula  approach  Is  utlLlzed;  Indications  of  how  to  proceed  in  the 
Lagranglan  approach  are  Included  at  the  end  of  the  section. 

Let  h'(x),  x  =  {x,y}  denote  the  total  still  water  depth,  and  let 

h’  =  h ( x )  -  5 (x)  (2.1) 


where  h(x)  Is  a  slowly  varying  depth  satisfying  the  mlld-slope  assumption. 


Unfortunately ,  no  laboratory  data  exists  to  test  these  hypotheses; 


rlflcation  of  the  present  results  and  conjectures  thus  require  further 
dependent  effort. 


.4  Coupled  Parabolic  Equations  for  Forward  and  Back-scattered  Waves 
We  now  consider  the  development  of  coupled  parabolic  equations  for 
rward  and  back-scattered  waves,  following  the  results  of  Radder  (1979)  and 
u  and  Tsay  (1983).  The  goal  is  to  obtain  an  extension  to  the  refraction 
suits  of  Mei  (1984)  (as  in  3.5  and  3.6)  to  cover  cases  where  6  varies 
•bitrarily  in  x  and  y,  and  where  y-variations  in  6  or  h  may  induce 
ifficiently  strong  amplitude  variations  to  warrant  the  introduction  of 
.ffraction  effects.  We  take  x  to  correspond  to  a  principal  propagation 
Lrection  and  assume  that  deviations  from  this  direction  are  small, 
jglecting  time  dependence  in  the  wave  amplitude,  (2.11)  may  be  written  in 
Lliptic  form  as 


p$  +  p  <j>  +  k  CC 

XX  X  X  g 


(CC  ( 6<j>  )  =  0 

g  y  y  k  Yy'y 


(4.1) 


lere 
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CC„ 


4ojW  '  6 
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+  0(k6)2  . 
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Figure  II. 5.  Reflection  from  a  patch  with  two  Fourier  components 

6  =  l)jslnAx  +  IHslnC  1 5Ax/8) ,  n=4. - ,  Dj/h  =  0.32, 

02 /h  =  0.0;  "  ,  D i /h  =  0.0,  D2/h  =  0.32; 

D ! /h  =  D2/h  =  0.226. 


For  case  one,  we  use  the  data  of  Figure  II. 3  and  set  m»2  so  that  the 
second  wave  has  half  the  wavelength  of  the  first.  Three  situations  are 
plotted  in  Figure  II. 4;  Dj/h  =  0.32,  D2/h  =  0.0  (as  in  Figure  II. 3);  D^/h  - 
0.0,  D2/h  =  0.32  (all  variance  at  the  shorter  wavelength);  and  D^/h  =  D2/I1  =* 
0.226  (variance  evenly  divided  between  each  component).  The  composite  bottom 
is  seen  to  produce  the  same  zeroes  in  R,  and  a  peak  in  R  associated  with 
resonance  with  each  component  of  the  bottom  is  apparent.  The  two  resonance 
peaks  remain  separate  and  clearly  distinguishable. 

For  the  second  case,  we  use  n*4  and  m=15/8,  so  that  the  second  component 
has  7  1/2  wavelengths  in  the  ripple  patch.  Curves  of  R  for  the  same  three 
distributions  of  ripple  amplitudes  as  described  above  are  given  in  Figure 
II. 5.  Now  the  reflection  coefficients  associated  with  each  component  acting 
separately  have  different  zeroes.  The  two  patterns  together  interact  and 
destroy  the  zeroes;  the  resulting  curve  of  R  varies  smoothly  over  the  range  of 
2k/i.  considered.  The  peaks  associated  with  resonant  scattering  from  each 
component  are  again  distinct  and  strong. 

The  strength  of  resonant  reflection  from  an  organized  barfield  implies 
that  a  broad  spectrum  containing  a  band  of  wavelengths  which  are  nearly  or 
exactly  resonant  with  the  bottom  may  experience  fairly  significant  reduction 
In  energy  density  near  the  resonance  wave  number.  Further,  the  reflected  wave 
field,  which  may  assume  the  form  of  a  fairly  narrow  spectral  band  around  the 
resonant  wave  number,  may  be  much  "groupier"  than  the  incident  wave  field. 

This  groupiness  may  lead  to  significant  forced  long  wave  motions  propagating 
in  the  offshore  direction.  These  motions  would  be  similar  in  form  to  the 
offshore  propagating  long  wave  caused  by  wave-group  pumping  of  the  surf  zone 
(studied  recently  by  Symonds  and  Bowen,  1984,  and  Symonds,  Huntley  and  Bowen, 
1982),  but  do  not  require  the  incident  wave  field  to  be  distinctly  groupy  in 


nature 
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2k/X 


Figure  11.3.  Reflection  coefficient  for  waves  normal 1  •  Incident  on  a 

sinusoidal  patch.  Case  2:  D/h  **  0.32,  n  *  4.  numerical 

results;  •  laboratory  data  from  Davies  and  Heathershaw  (1984) 
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Reflection  coefficient  for  waves  normaily  incident  on  a 
sinusoidal  patch.  Case  1:  D/h  =*  0.16,  n  *  10.  numerical 

results;  •  laboratory  data  from  Davies  and  Heathershaw  (1984) 


are  given  In  Figures  11.2  and  11.3  for  cases  1  and  2.  respectively  In 
comparison  to  the  laboratory  data.  Solutions  were  obtained  using  a  grid 
spacing  of  Ax  -  1/20  In  order  to  obtain  accurate  results  at  large  2k/X.  The 
numerical  results  are  nearly  indistinguishable  from  the  analytic  results  of 
Davies  and  Heathershaw  (1984)  with  the  exception  of  an  occasional  downshift  In 
the  positions  of  the  peaks  and  zeroes  of  R  by  an  amount  2k/X  ~  0.01.  This 

shift  may  be  due  to  the  neglect  of  nonpropagating  waves  at  the  ripple  patch 
edges  as  well  as  the  numerical  errors  involved  In  the  computation  scheme.  The 
results  are  based  on  the  same  Information  as  In  Mel's  approach,  and  compare 
well  to  his  predictions  over  the  range  of  validity  of  the  theory  for  near¬ 
resonance. 

Laboratory  data  from  Davies  and  Heathershaw  are  included  In  Figures  II. 2 
and  II. 3  for  comparison  to  the  numerical  results.  Davies  and  Heathershaw 
provide  Indications  of  how  much  the  data  at  low  R  Is  contaminated  by 
reflections  at  the  end  of  the  wave  channel;  no  effort  has  been  made  to  correct 
the  data  for  this  effect. 

As  a  second  example,  we  consider  the  reflection  of  waves  from  a  bed 
formed  by  the  superposition  of  two  Fourier  components. 

6  =  DjslnXx  +  D2slnmXx  0  x  <_  2wn/X 

We  consider  two  cases;  one  where  both  Fourier  components  have  an  Integral 
number  of  wavelengths  In  the  ripple  patch,  and  the  other  where  the  ripple 
patch  terminates  at  a  half  wavelength  for  the  second  component.  These  cases 
differ  since,  for  case  one,  the  zeroes  of  the  reflection  coefficient  for  each 
component  separately  coincide;  while  for  the  second  case  the  zeroes  for  each 
component  do  not  coincide. 


6  -  Dsin(Ax)  0  <  x  <  nl 


(3.11) 


where  Z  ■  2n/A  ts  the  bed  wavelength  and  n  la  the  number  of  bed  ripples.  The 
Bragg-scatterlng  condition  corresponds  to  2W/A  -  1.  Equation  (3.8)  Is  written 
In  finite  difference  form,  and  radiating  boundary  conditions  are  applied 
according  to 


>x  *  ~  ik($  “  2#x)  Xj  <  0 


(3.12) 


♦  *  lk$ 

x 


x2  >  nl 


(3.13) 


where 


*1  =  G 


(3.14) 


Is  the  Incident  wave  of  unit  amplitude.  The  resulting  tridiagonal  matrix  Is 
Inverted  using  a  double-sweep  algorithm. 

We  first  consider  the  experimental  results  of  Davies  and  Heathershaw  and 
compare  numerical  results  to  the  cases 


1.  n  -  10  ,  D/h  =  0.  16 

2.  n  =  4  ,  D/h  *  0.32 


Computed  resuLts  for  reflection  coefficient  R  corresponding  to  a  reflected 


—  i  lev 

$K  “  Re  ;  x  <  0 


(3.15) 


/ '.'.v  V 
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Met.  All  of  Mel's  subsequent  results  nay  be  obtained  f roe  (2.11)  following 


this  procedure  and  his  assuaptions. 

For  later  use,  we  define  an  unsealed  f requency-llke  tera  Q'  according  to 

no  -  ft’(kD)  (3.7) 

where  ft'  ~  0  (1).  ft'  will  be  advantageous  in  the  general  case  since  it  Is  not 
dependent  on  the  geometry  of  the  undulations.  (Note  that  ft*  aay  still  be  a 
slowly  varying  function  of  the  mean  depth.) 

II. 3. 2  One-Dtmens tonal  Reflection  from  a  Ripple  Patch 

We  now  test  the  reduced,  elliptic  form  of  (2.11)  for  the  case  of  waves 
normally  Incident  on  a  finite  ripple  patch.  Variations  in  the  y-dlrectlon  are 
neglected;  the  resulting  problem  Is  equivalent  to  that  studied  by  Davies  and 
Heathershaw  if  we  restrict  attention  to  sinusoidal  undulations  of  constant 
amplitude  and  constant  mean  depth  h.  After  setting  h  »  constant,  (2.11) 
reduces  to 

♦  xx  +  ^2d»  -  ')(,s,i>x)x  “  0  (3-8) 

g 


where 

~ 

♦  =  <t>(x)  e  ,  (3.9) 

ft'  Is  defined  by 

ft'  -  - y  '  (3.  10) 

4oi  cosh  kh 

following  (3.7),  and  where  we  have  negLected  time  dependence  In  the  wave 
amplitude.  For  the  case  of  sinusoidal  bed  oscillations,  we  may  take 


We  consider  an  incident  wave 


“  2?  A(x>t)  e 


l(kx-wt) 


and  reflected  wave 


♦  •  -  U  B(x*t)  e 


i(-kx-wt) 


0.1) 


(3.2) 


where  A  and  B  are  slowly  varying  functions  of  x  and  t.  The  conditions  for 
resonant  Bragg-  scattering  are  satisfied  when  the  bottom  undulation  has  a 
wavelength  of  one-half  the  surface  wavelength,  or 

X  =*  2k  (3.3) 

Employing  this  assumption,  we  proceed  by  assuming  that  derivatives  of  A  and  B 
are  0(k6)  in  comparison  to  A  and  B  and  keep  terms  only  to  0(k6).  Defining 

-  - Sk-°£-  (3.4) 

cosh  xh 

1  Wx 

and  collecting  terms  of  like  powers  in  e  ,  we  obtain 


A  +  C  A 
t  g  x 


-  W0B 


(3.5) 


B 


t 


CB 

(?  x 


-  m  a 

o 


(3.6) 


which  Is  a  special  case  of  the  results  of  Mel,  neglecting  mean  bottom  slope 
and  oblique  angle  of  incidence.  Solutions  of  (3.5)  and  (3.6)  In  relation  to 
the  experimental  results  of  Davies  and  Meathershaw  are  discussed  in  detail  by 


1*  ■  /  Ut  ♦  Y  (7h9)2  ♦  2"  <*z>2  *  8Z1  dt 

and  Chen  expanding  Che  InCegral  about  z  ■  -  h.  Equations  (2.10)  and  (2.11) 

2 

follow  after  collecting  teras  of  0(ka)  and  partially  Integrating  the 
variation  of  the  integral  of  L  over  the  propagation  space,  as  In  Kirby  (1984). 

II. 3  Correspondence  to  Previous  Results 

In  order  to  demonstrate  the  completeness  and  generality  of  the  wave 
equation  (2.11),  we  first  consider  a  reduction  of  the  equation  to  the  coupled 
evolution  equations  of  Mel  (1984)  for  resonant  Bragg-acattering  by  a  finite 
patch  of  ripples.  We  then  employ  a  finite  difference  form  of  (2.11),  after 
neglecting  time  dependence,  to  study  the  one-dlmenslonal  (x  only)  reflection 
for  a  range  of  incident  wavelengths,  and  compare  our  results  to  the  data 
presented  by  Davies  and  Heathershaw  (1984). 


II. 3.1  Resonant  Bragg-scatterlng 

Consider  the  particular  example,  studied  by  Davies  and  Heathershaw  (1984) 
and  Mel  (1984),  of  waves  propagating  normally  over  a  ripple  patch  extending 
uniformly  to  +_  ™  in  y.  Depth  h  is  taken  to  be  constant,  while  6  Is  given  by 


6 


-iAx. 
e  ) 


0  <  x  <  L  , 


where  0  £  x  <  L  Is  the  range  of  the  ripple  patch,  and 


6 


The  integrals  are  manipulated  to  finally  obtain 


*tt~  vh*(cVhJ)  *  (“*2  -  k2cV> ♦  *f2  ihV(6V> 
0  , 

+  g  /  fV^f  dz*  ♦  gfVhf*7hh  I  * 

-h  “h 

+  gfVhf.Vh6  |  *  +  g6(fV*f)  |  « 

-h  “h 


+  2g6fVhf«Vh#  |  -  0  (2.10) 

- h 

where  C  *  w/k  and  C~  -  3ui/3k.  The  last  five  terms  are  either  proportional  to 

O 

(V^h)^  or  <5Vhh,  and  are  thus  second  order  In  the  small  parameter.  Neglecting 
them  and  substituting  for  f(z  -  -  h)  then  gives 


-  V  .(CC  7,*)  +  (w2  -  k2cc  )?  + 

h  g  h  g 


cosh  kh 


V<4V> 


«  0(k5)2 


(2.11) 


Equation  (2.11)  governs  the  value  of  the  potential  at  the  free  surface  for  an 
arbitrary  wave  motion.  Neglecting  the  term  in  6  yields  a  time  dependent  form 
of  Berkhoff's  (1972)  equation  for  the  slowly  varying  bottom  alone. 

This  completes  the  derivation  using  the  Green's  identity  method.  The 
Lagrangian  formulation  proceeds  by  writing  the  integral  of  the  total  pressure 


over  the  local  depth 


Vhh/kh  «  1  ,  -  { 3/3x  ,  3/3y } 


(2.2) 


6(x)  represents  rapid  undulations  of  the  depth  about  the  mean  level,  as 
Indicated  In  Figure  II. 1.  We  consider  the  problem  to  be  linearized  In  wave 
amplitude  but  retain  first  order  terms  In  the  bed  undulation  amplitude.  We 


0(Vhh/kh)  ~  0(k6)  «  1 


(2.3) 


Linearizing  the  free  surface  boundary  conditions  and  expanding  the  boiton 
boundary  condition  about  z  -  -  h,  we  obtain  to  0(k6) 


V  +  *zz  =  0  5  -  h  <  2  1  0 


(2.4) 


*tt  +  glpz  *  0  }  2  “  0 


(2.5) 


~  7hh*Vh^  +  Vh’^67h^  ’  2  *  ~  h 


(2.6) 


To  leading  order,  the  solution  to  (2.4)  -  (2.6)  may  be  expressed  as 


$(x,z,t)  ’  f(x,z)  <Kx.O  +  E  non-propagating  modes 


(2.7) 


where  f  =  cosh  k(h+z)/cosh  kh  Is  a  slowly  varying  function  of  x,  and  where 


w  =  gk  tanh  kh 


(2.8) 


locally,  with  ^  being  the  fixed  angular  frequency  and  k  the  wavenumber.  We 


then  use  Green’s  second  identity  to  extract  the  propagating  component  of  4> ; 


where,  Co  leading  order  In  (k5),  p  *  is  given  by 

p~l  =  (CCg)-1  { l  +  4  (^— )  6+  0(k5)2}  (4.5) 

g 
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Next  we  denote  an  operator  y  $  according  to 

y2*  =  k2CC  p~S  +  p^CC  *  )  -  (6*  ) 

g  g  y  y  k  y  y 

or 


Y2*  =  k2  {(l  +  4  (—■)  6)  *  + 


2 

k4CC 


(CC  ♦  )  -  ^ 

g  y  y 


(r)  (6S)y» 


+  0(k6)2 


(4.6) 


The  corresponding  pseudo-operator  y4>  Is  obtained  by  the  expanding  the  square 
root  to  give 


Y$=k{(l+2  )  «)♦  + 


2k2CC 


(CC  $  ) 

g  y  y 


kT  (I1)  <  Vyi 

k  g 


+  0(k5  y 


(4.7) 


We  now  follow  a  simple  scheme  for  obtaining  the  coupled  parabolic 
equations.  Let  $  be  expressed  as  the  sum  of  forward  and  backward  propagating 


waves ; 


We  then  assume  the  coupled  equations: 


“  iY$  +  F($  ) 


-  -m  -  f($+,$~) 


(4.9) 


where  the  coupling  term  F  Is  unknown.  Repeated  substitution  of  (4.9)  In  (4.4) 


finally  yields 


F(*+,*  ) 


(YP)x 


(4>  -  4>  ) 


(4.10) 


which  may  be  expanded  to  give 


(kCC  ) 

'(♦*•♦')  --(^cP-(r)  s,t 

g  g 


(4.11) 


to  leading  order  In  (k6).  The  coupled  equations  In  expanded  form  are  given  by 


»*  -  Ik  (1  ♦  2  (§4  6)  ♦*  .  2kCC~  <CC,*y)y 


01  o.  ^  (kCC  ) 

-  i r  £-)  K>y  -  -lEcP  -  ♦’> 

g  g 


+  (I1)  «x(/-  <0  )  ; 

g 


(4. 12) 


♦;  =  -  ik  (i  +  2  Cr)  «} 


2^  <cWy 


9«  o.  (kCC  > 

+  r  (r)  <  Vy  +  ~2kCC~ <♦  -  ♦’) 
g  g 


-  (jf-J  «x(*f-  ♦’)  • 

g 


(4.13) 


a*  L*  %-*  O  ^  a.*  a.'  •  **  mJ*  •*. 


We  introduce  the  complex  amplitudes  A,B  according  to 


i  Ik  x 
*+  =  -  ^  Ae  ° 

U) 


,  -ik  x 
_  iiL  -  o 


4>  - - “■  Be 

w 


where  kQ  is  a  reference  wavenumber.  (4.12)  and  (4.13)  become 


2ikCC  A  +  {2k  (k-k  )  CC  +  i(kCC  ) 
g  x  o  g  g  x 


+  2u2 '  [ 2k6  -  16  ]}  A  +  (CC  A  )  -  — ~  (6A  ) 

x  1  g  y  y  k  y  y 

-2ik  x 

=  { i (kCC  )  -  2iwJl'6  }  Be  °  (4.14) 

g  x  x 


2ikCC  B  +  {-  2k (k-k  )  CC  +  i(kCC  ) 
g  x  1  o  g  g'x 


-2wtt'  [ 2k6  +  i6  ]}  B  -  (CC  B  )  +  (6B  ) 

xJ)  g  y  y  k  y  y 


2ik  x 

=  { i(kCC  )  -  2iwft'5  }  Ae  °  (4.13) 

^  X  X 

The  correspondence  between  the  present  equations  and  the  results  in 
Section  3.1  may  be  seen  by  substituting  for  6  using  (3.11)  and  neglecting 
terms  with  rapidly  oscillating  coefficients.  Note  that  the  equations 
developed  by  Mei  neglect  diffraction  effects  as  well  as  the  coupling  between 
the  forward  and  backward  propagating  waves  over  the  slowly  varying  depth.  T 
present  equations  Include  these  effects.  Further,  they  reduce  to  a  set  of 
equations  equivalent  to  those  of  Liu  and  Tsay  (1983)  when  6  Is  neglected. 


In  the  following  example,  (4.14)  -  (4.15)  are  discretized  according  to 
the  Crank-Nicolson  method.  The  solution  technique  is  equivalent  to  that  used 
by  Liu  and  Tsay;  hence,  the  details  are  omitted  here. 

We  construct  a  two-dimensional  patch  of  ripples  of  finite  extent  in  the  x 
and  y  directions.  Ripples  with  length  i  are  alligned  with  crests  parallel  to 
the  y-axis.  The  patch  is  symmetric  about  the  x-axis  and  has  overall 
dimensions  of  nil  and  2n l  in  x  and  y,  where  n  is  the  number  of  ripple 
wavelengths.  The  topography  is  given  by  h  *  constant  and 

f Dsin(Xx)cos(Xy/4n) ; 

<5(x,y)  =  <  |x|<n£/2,  |y|<n£  (4.16) 

to  ;  |  x  |  >n£ /2,  |yj>n£ 

The  computational  domain  is  given  by  -3  <  x/£  <  3  and  0  <  y/l  <•  6.  We 
consider  ripples  similar  to  those  of  Figure  II. 3,  with  D/h  =  0.32  and  n=4. 
Results  were  computed  for  the  resonant  case  2k/X  *  l  and  are  plotted  along 
with  the  bottom  contours  in  Figures  II. 6a  (for  incident  amplitude  | A | )  II. 6b 
(for  reflected  amplitude  |b|),  and  Figure  II. 6c  (for  the  total  wave  field). 

The  present  results  were  obtained  using  four  Iterations  of  the  forward- 
backward  calculation,  which  was  sufficient  to  provide  a  reasonable  degree  of 
convergence.  As  in  Section  3,  a  grid  spacing  of  Ax  =  Ay  =  X/20  was  used. 
Unfortunately,  no  laboratory  data  exists  to  test  the  two-dimensional  model; 
verification  was  limited  to  checking  that  results  of  the  coupled  parabolic 
model  are  equivalent  to  results  using  the  elliptic  model  for  the  one- 
dimensional  cases  studied  in  Section  3. 


a)  Incident  wave  field  | A/ AQ | 

Amplitude  contours  with  respect  to  Incident  wave  amplitude 
waves  propagating  over  two-dimensional  ripple  patch. 

-  bottom  contours;  ■  -  amplitude  contours 


*  l  ■  ^ 


II. 5  Waves  Trapped  Over  Long  Ripple  Patches 

As  a  final  example,  we  consider  the  propagation  of  a  trapped  wave  along  a 
ripple  patch  which  extends  to  y  »  t  *.  This  case  has  possible  applications  to 
the  nearshore  zone,  where  waves  may  be  directed  onto  a  nearshore  bar  system  by 
reflection  from  headlands  or  coastal  structures.  This  problem  has  been 
investigated  previously  in  the  context  of  the  mild-slope  approximation  by 
Lozano  (1977);  the  present  analysis  is  presumably  applicable  to  steeper  bar 
systems  than  would  be  allowed  for  under  the  mild-slope  assumption.  Further, 
in  the  case  of  long  coastally-trapped  waves,  the  bar  system  may  serve  as  a 
wave  guide  for  trapped  wave  energy;  Kirby  _£t__al_  (1981)  have  suggested  that 
this  trapping  may  effectively  stabilize  the  offshore  distribution  of  nodes  and 
antinodes  of  infragravity  edge  waves,  rendering  the  bar  system  less  subject  to 
destructive  effects  of  a  slowly  changing  wave  climate. 

The  general  case  of  trapped  wave  motion  over  a  variety  of  offshore 
profiles  may  be  handled  in  a  manner  similar  to  that  provided  by  Kirby  et  al 
(1981);  we  remark  that  this  method  has  also  been  used  recently  by  Mclver  and 
Evans  (1984)  to  study  the  waves  trapped  above  a  submerged  horizontal  cylinder 
located  close  to  the  water  surface. 

We  denote  l  as  the  longshore  (y)  component  of  the  wave  number  vector  and 
assume  that  h  and  6  are  functions  of  x  only  (no  longshore  variation  in  the 
depth  profile).  Letting 

♦  (x,y)  =  iKx)eliy ,  (5.  1) 


we  obtain  the  second  order  ordinary  differential  equation  governing  i|>(x); 

4ojn'<5 


(cc  -  6-}  \>  +  {(cc )  - 

1  g  k  J  xx  1  g  x 


ir^x 


M(X2  -  i2)ccg  +  =  0 


(5.2) 
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This  then  provides  an  eigenvalue  problem  for  i|>  and  l  which  must  be 
solved  together  with  suitable  boundary  conditions.  As  a  first  example,  we 
consider  the  simple  case  of  a  rectangular  ridge  of  height  D  and  width  2L 
resting  on  a  flat  bottom  of  depth  hj.  The  depth  of  water  over  the  ridge  is 
given  by 


h2  -  hL  0  <  h[  (5.3) 

The  geometry  is  shown  in  Figure  II. 7.  This  case  is  of  interest  because  it  may 
be  solved  analytically  both  by  the  present  method  as  well  as  by  using  the  mild 
slope  equation  in  its  basic  form  (see  Smith  and  Sprinks,  1975,  for  the 
relevant  formulation  and  matching  conditions).  In  addition,  the  shallow  water 
limits  of  both  solutions  may  be  compared  to  the  results  presented  by  Mei 
(1983,  Section  4.6).  This  case  then  provides  us  with  a  means  for  evaluating 
the  approximate  upper  limit  on  D/h^,  which  is  determined  by  checking  for 
divergence  of  the  solution  for  the  present  theory  from  the  mild-slope  theory 
as  D/hj  becomes  large. 

For  the  case  of  constant  depth  h2  =  h1  -  D  over  the  ridge  and  hj  away 
from  the  ridge,  we  divide  the  fluid  domain  into  three  regions  and  obtain 

*(2)  +  (k2  -  l1)  i>(2)  =0  ;  x  >  L  (5.4a) 

XX  1 

■4-(3)  +  (k2  -  l2)  <p(3)  =0  ;  X  <  -  L  (5.4b) 

xx  1 

and 

(k2  -  aH2)  4’^^  =  0  ;  -  L  <  x  <  L 

where  kj  is  the  wavenumber  corresponding  to  h^  and 


(5.4c) 


a  =  1  "  4  Cc”)  D 

g 


(5.5a) 


=  1 


2kiN 


n  sinh2k^ 


(5.5b) 


In  order  for  trapped  waves  to  exist,  solutions  must  be  exponentially  decaying 

2  2 

away  from  the  ridge  in  regions  2  and  3.  This  requires  that  l  >  k^  while  a 
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separate  condition  on  '  implies  that  l  <  k^/a.  The  solutions  for  each 
region  are  given  by 


,j/2>  =  Ae"  ra(x-L) 

(5.6a) 

4,0)  =  Bem(x+L) 

(5.6b) 

*(!>  =  CeUx  +  De_1Xx 

(5.6c) 

where  A,  B,  C,  D  are  arbitrary  constants, 

1/2 

m  =  (JT  -  kp  >0  ,  (5.7) 


and 


X 


(k2/a  -  l2) 


1/2 


>  0 


(5.8) 


Matching  conditions  are  provided  by  considering  the  Integral  of  the  governing 
equation  across  the  jump  discontinuity  at  the  ridge  boundary;  we  get 


x-L 


(5. 9a) 


x--L 


,<l>  (2) 

a\p  =  si; 
x  x 


x  =  L 


(5.9b) 


(5.9c) 


,(1)  ,(3) 

ctU>  =  i b 
x  x 


x  =  -  L 


(5.9d) 


The  solution  to  (5.6  -  5.9)  may  be  found  in  a  straightforward  manner;  the 
notation  of  Mei  (1983)  is  used  to  facilitate  comparison  to  the  shallow  water 
case.  Trapped  waves  are  given  by  the  solutions  (for  £ )  of  the  relation 


U*2  -  £2) 


cot£  n  even 


-tan£  n  odd 


(5.10) 


where  n  is  the  mode  number  of  the  trapped  wave.  Here,  £  and  £*  are  given  by 


£  =  {(kj/a  -  l2)  1  L 


(5.11) 


£*  =  (k  1L)  ((1  -  a)/a) 


(5.12) 


The  condition  that  £  <  £*  from  (5.10)  provides  an  upper  cutoff  limit  for  n;  n 
Increases  for  increasing  L  or  decreasing  (increasing  D). 

Following  Smith  and  Sprinks  (1975),  the  mild-slope  equation  may  be  solved 
in  analogous  fashion.  The  solutions  are  given  by  (5.10)  with  the  newly 
defined  coefficients 


v  \ 


(5. 13a) 


a  =  (CCg)2/(CCg)l 

„  „  1/2 

?  2 

x  =  (kj  -  *  ) 


5  =  {(k2  -  Z2)  }  L 

1/2 

5*  =  {(k2  -  k2)  }  L 


(5. 13b) 

(5.13c) 
(5. 13d) 


or  the  case  of  the  present  geometry,  where  the  mild-slope  equation  is  valid 
xcept  across  the  jump  discontinuities,  we  may  regard  (5.10)  together  with 
5.13)  as  a  "correct"  solution  and  compare  those  results  to  the  present,  smal 
>ottom  amplitude  solution  given  by  (5.10)  together  with  (5.8  -  5.12). 
Comparing  the  various  coefficients,  we  see  that  the  present  theory  should 

ipproximate  the  general  depth,  mild-slope  theory  well  in  regions  where  a  from 

1/2 

(5.5)  is  a  good  approximation  to  a  from  (5.13a),  and  where  k^/a  '  (a  from 

(5.5) )  is  fairly  close  to  k2. 

1/2 

Plots  of  the  two  a's  and  of  the  ratio  (kj/ct  '  )/k2  are  given  in  Figures 
11.8a  and  b  for  a  range  of  t^/hj  and  kjhj  values.  Interestingly,  the 
deviations  between  the  two  theories  are  greatest  in  the  intermediate  depth 
range  and  apparently  decrease  in  deepwater.  This  result  is  due  to  the  fact 
that  the  depth  over  the  step  becomes  effectively  deep  unless  h2/hj  is  small, 
rendering  the  step  ineffective  in  trapping  or  reflecting  waves.  In  shallow 
rfater,  both  theories  give  solutions  which  are  asymptotic  to  the  results  given 
by  Met  (1983).  The  asymptotic  values  of  a  and  £*  are 


a  -  h2/hL 


1/2 


L(l“h2/hl) 


1/2 


h  *  0 


(5.14) 


The  two  theories  are  essentially  different  in  intermediate  depths; 
however,  the  results  differ  little  over  the  range  0.6  <  h~/h,  <  1.0,  or, 


igure  I [.8.  Comparison  of  parameters  for  trapped  waves  over  a  rectangular 
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:ter  eliminating  n j  according  to  (2.9).  Variation  of  L ^  with  respect 
>  4>.^  and  partial  integration  yields  a  continuity  equation 


*2  *  Vt  +  Vh 


(hvt-)  --vh 


(l>  V  4  )  +  V 
^1  hTr  h 


-  0 


(2.14) 


■\  ^  +  b0)  may  be  eliminated  between  (2.13)  and  (2.14)  to  yield  a  forced  wave 
quation  for  the  quasi-steady  motion  after  averaging  over  the  phase. 

The  equation  containing  the  nonlinear  modification  to  the  linear  wave 
quation  (2.11)  may  be  obtained  by  varying  +  L^  together  by  n  j  and  4>  and 
liminat tng  np  yielding 


e  {$!  -  V(CCgVl)  +  (u>2  "  k2CC  )$i}  =e3(N.L.T.}  (2.15) 

N.L.T.}  is  a  complicated  expression  involving  products  of  c*>  ^ ,  (q  2  +  t>2 ) » 

2  and  >2»  and  Is  given  in  Appendix  IIL.C  for  completeness.  We  remark  that 
erms  in  (N.L.T.}  have  been  manipulated  by  substitution  using  the  linearized 
elat ionships ,  in  analogous  form  to  the  treatment  of  higher  order  terms  in  the 
oussinesq  wave  formulation.  Equation  (2.15)  is  in  the  form  of  a  second  order 
yperbollc  equation  for  a  general  wave  motion  in  {x,y,t};  however,  {N.L.T.} 
ontains  components  proportional  to  the  thLrd  harmonic  These  terms  are 

liminated  in  the  derivation  of  the  evolution  equations.  An  equation  similar 
o  (2.15)  has  been  given  previously  by  Liu  and  Tsav  (1984)  for  the  case  of 
rogresstve  waves. 
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at  parameter,  It  Is  clear  that  the  linear  equation  governing  the  behavior  of 
parameter  of  0(en)  will  come  from  the  variation  of  with  respect  to  that 


iraraeter.  Thus  the  linear  equation  governing  <|>  ^  will  be  contributed  by  L^, 
llle  L4  will  contribute  the  0(e  )  terms  In  phase  with  4^. 

We  now  derive  governing  equations  which  are  generally  applicable  to  any 
Dtlon  described  by  the  potential  (2.4a)  and  contribute  nothing  in  this 
pproxlmatlon.  The  linear  wave  equation  is  derived  by  first  varying  L2  with 
espect  to  H[,  yielding 


n 


1 


(2.9) 


he  free  surface  boundary  condition.  Varying  L2  with  respect  to  $  ^  and 
•erforming  a  partial  Integration  over  the  propagation  space  yields 


<Ii,iVi>+  vV0  (2*I0) 

Eliminating  Hj  between  (2.9)  and  (2.10)  and  inserting  the  values  of  the 
integrals  yields 

4>.  -V  •  (CC  V  I)  +  (w2  -  k2CC  )♦.  =  0  (2.11) 

o  **  1  ri 

»hlch  Is  the  time  dependent  form  of  the  mild  slope  equation  of  Berkhoff 
( L 97 2 ) .  Here, 

2 

C  =  w/k  ;  Cg  =  3w/9k  ;  u  =gktanhkh.  (2.12) 


Varying  with  respect  to  (n  2  +  bj)  yields  a  free  surface  boundary 
:ondltlon,  which  may  be  written  as 
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L=|-(n"+h)+el1<t>1  +  e  { 1 2<J)  2  +  (h+n)(4>2  -  Y  2>} 
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(V2> 
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2  z  (V  3  z  ~  ~  4  z 

+  *  II,  1  "1  '  1 1 , 2  *1*2  +  £  X2,2  ~2 


(2.6) 


The  integrals  I  are  over  the  total  depth  and  are  defined  in  Appendix  III. A. 
The  I’s  may  be  expanded  about  z=0  in  Taylor  series  according  to 


I  =  I'  +  en  LI  ’ '  +  e2{  (n2+b2)I"  +  +  e  3{  2n  2+b2)I "  • 


+nL3IIV}+  0(P)  (2.7) 


as  shown  in  Appendix  III. A.  Substituting  the  expansions  In  (2.6),  expanding 
the  remaining  appearances  of  n,  and  retaining  terms  to  0(eP)  leads  to  the 
expression 


L  =  L0  +  eLj^  +  e2L2  +  e3L3  +  e4L4 


(2.8) 


The  individual  coefficients  are  given  in  Appendix  III.B.  This  ordering  has 
also  been  utilized  by  Dysthe  (1974)  after  averaging  L  over  the  phase  function 
according  to  Whltham's  method.  Serveral  properties  of  (2.8)  with  respect  to 
the  desired  governing  equations  can  be  mentioned.  Since  variation  with 


respect  to  a  dependent  parameter  will  reduce  the  order  of  L  by  the  order  of 


6  /  /  L  dx  dt  =  0 


(2.3) 


where  the  Integration  is  over  the  propagation  space  {x,y,t}. 

Since  we  wish  to  study  the  spatial  evolution  of  a  time-independent  wave 

field,  it  is  sufficient  to  choose  forms  for  4>  and  h  based  on  a  s low-modulation 

solution  of  the  governing  equations,  as  in  Whitham  (1967).  Consequently,  we 

choose  a  wave  steepness  scale  e  and  modulation  scale  u  and  propose  a  priori 
2 

that  w~e  •  This  assumption  produces  a  mild  slope  approximation  in  which 
bottom  slope  terms  and  nonlinear  terms  are  isolated  from  each  other,  as  in 
Kirby  and  Dalrymple  (1983).  We  choose  4>  and  n  according  to 


$  =  ef  1(ux12H1(x,t)  +  e  {f 2(ux,z)<t>2(x,t)  +  4> 2  (x.^ t)-Jy 2 (v*2L» t)dt} 


2 

n  =  en1  (x,t)  +  e  (n2(x,t)  +  b2(x,pt)} 


(2.4a) 

(2.4b) 


Here,  4>  1  >9  j  »$2 ,r|  2  corresPon<^  to  first  and  second  order  wavelike  components, 

<)>2  is  the  potential  for  wave-induced  mean  flow,  b£  represents  the  wave-induced 
set-down,  and  y 2  is  related  to  the  Bernoulli  constant.  We  remark  that  in  a 
partial  standing  wave  Y2  cannot  be  trivially  eliminated  by  choice  of  b2,  and 
that  ^2*  b2  and  ^2  ma-*r  ^ave  ^ast  variations  as  indicated.  The  quantities  fj 
and  f2  ate  given  by 


coshk(h+z) 
3  — - •  “  ■ 

1  coshkh 


cosh2k(h+z ) 
sinh^kh 


h  =  h(wx) 


(2.5) 


Substituting  (2.4)  in  (2.1)  and  performing  the  integration  results  in  the 
expression 
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III. 2  Derivation  of  the  Equations  Governing  Wave  Propagation 

We  wish  to  derive  the  equations  governing  the  forward  and  backward  moving 
wave  components  in  a  partial  standing  wave,  where  x  is  taken  as  the  (positive) 
direction  of  travel.  Previous  derivations  of  the  nonlinear  Schrodinger 
equation  for  forward  scattered  waves  alone  have  relied  on  the  WKB  formulation 
and  a  multiple  scale  expansion  of  the  governing  equations.  However,  a  scheme 
for  providing  coupled  equations  has  not  been  devised  using  that  approach. 

Here,  we  will  rely  instead  on  a  variational  formulation  using  the  Lagrangian 
for  irrotational  motion  of  an  inviscid  fluid,  given  by  Luke  (1967).  After 
deriving  equations  for  a  general  wave  motion  in  two  horizontal  dimensions 
(x,y)  and  time,  the  results  will  be  specialized  to  the  case  of  two  waves 
propagating  in  an  anti-parallel  direction.  The  method  used  here  is  further 
discussed  in  Kirby  (1983)  and  Kirby  and  Dalrymple  (1984)  in  connection  with 
the  problem  of  wave-current  interaction. 


III. 2.1  The  Lagrangian  Formulation  and  Governing  Equations 

The  Lagrangian  for  irrotational  motion  is  given  by  (Luke,  1967) 


L  =  /  p { j  (vh+)2  +  J  ('JO  2  +  dz 
-h 


(2.  1) 


Here,  n(x,y,t)  is  the  instantaneous  position  of  the  water  surface  with  respect 
to  still  water  level  z=0,  and  h(x,y)  is  the  local  water  depth.  The  potential 
i>(x,y,z,t)  is  related  to  the  fluid  velocity  according  to 


(2.2) 


,3  3  , 

Further,  7.  denotes  a  horizontal  gradient  vector  {  - —  ,  - —  }  and  subscripts 


denote  differentiation.  The  corresponding  variational  principle  is  given  by 


The  results  mentioned  above  are  based  on  the  assumption  that  the 
reflected  wave  Is  absent  or  negligible.  This  assumption  is  certainly  valid 
locally  In  a  slowly  varying  domain;  however,  a  sizeable  reflected  component 
may  accumulate  for  waves  propagating  over  long  distances  or  over  fairly  abrupt 
obstacles,  such  as  those  studied  in  Chapter  II.  Liu  and  Tsay  (1983)  have 
developed  an  iterative  scheme  based  on  coupled  equations  similar  to  those  of 
Radder  (1979),  and  have  shown  that  the  coupled  method  for  forward-  and  back- 
scattered  waves  is  capable  of  producing  results  in  agreement  with  a  finite 
element  solution  of  Berkhoff's  equation,  where  the  entire  wavefield  is 
calculated  simultaneously  (Tsay  and  Liu  (1982)). 

In  this  study,  we  extend  the  results  of  Kirby  and  Dalrymple  (1983)  and 
Liu  and  Tsay  (1983)  to  study  the  gradual  reflection  of  Stokes  waves  by  slow 
depth  variations.  In  addition,  the  results  of  Chapter  II  are  incorporated  in 
order  to  allow  for  the  application  of  the  model  to  reflection  from  sand  bars 
and  related  features.  In  Section  2,  a  derivation  of  the  equations  governing 
the  evolution  of  a  slowly  varying  partial  standing  wave  train  are  derived  in 
order  to  obtain  the  nonlinear  coupling  coefficients  between  incident  and 
reflected  waves.  Section  3  then  presents  a  modification  of  the  splitting 
approach  used  in  Chapter  II,  which  produces  the  needed  set  of  coupled 
parabolic  equations.  In  Section  4,  we  study  a  special  case  of  waves  normally 
incident  on  topography  varying  in  one  direction,  in  order  to  evaluate  the 
effect  of  neglecting  interaction  with  wave-induced  currents  which  enter  the 
equations  at  third  order.  The  parabolic  approximation  is  then  applied  to  the 
study  of  a  two-dimensional  problem  in  Section  5. 


Chapter  III 


On  the  Gradual  Reflection  of  Weakly-Nonlinear  Stokes 
Waves  in  Regions  with  Varying  Topography 


III. I  Introduction 

Various  recent  studies  have  explored  the  application  of  the  parabolic 
equation  method  to  the  problems  of  wave  propagation  in  a  slowly  varying 
domain.  In  the  linear  wave  approximation,  Berkhoff  (1972)  has  provided  a 
second  order  wave  equation  governing  the  surface  potential  <j>  (x,y)  of  a 
harmonic  wave,  based  on  the  assumption  of  modulation  length  scales  smaller 
than  the  wave  steepness  scale.  Radder  (1979)  obtained  coupled  parabolic 
equations  for  forward-  and  back-scattered  waves,  where  it  is  assumed  that  both 
waves  travel  at  a  small  angle  to  a  prespecified  direction.  A  model  for  the 
forward-scattered  wave  field  alone  is  then  obtained  by  neglecting  the  coupling 
terms.  Berkhoff,  Booij  and  Radder  (1982)  have  provided  a  data  set  for  wave 
amplitude  in  an  area  of  waves  focussed  by  a  submerged  shoal,  and  have  used  the 
data  to  test  Radder's  model.  Additional  results  for  forward-scattering  alone 
in  a  linear  approximation  have  been  provided  by  Liu  and  Met  (1976)  and  Tsay 
and  Liu  (1982).  Kirby  and  Dalrymple  (1983)  have  extended  the  linear  wave 
approximation  to  include  the  effect  of  cubic  nonlinearity  for  Stokes  waves, 
which  leads  to  a  nonlinear  Schrodinger  equation  for  the  wave  amplitude.  Kirby 
and  Dalrymple  (1984)  have  tested  the  nonlinear  model  in  comparison  to  the  data 
set  of  Berkhoff,  Booij  and  Radder  and  have  shown  a  distinct  improvement  in 
agreement  between  model  and  data  in  comparison  to  the  results  of  Radder's 
model.  The  nonlinear  model  has  thus  been  shown  to  be  a  relevant  addition  to 
the  study  of  combined  refraction-diffraction. 
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marked  contrast  to  the  case  of  edge  waves  on  a  plane  beach,  where  the 
shoreline  singularity  dominates  wave  amplification. 

II. 6.  Conclusions 

This  chapter  has  provided  a  general  wave  equation  for  linear  surface 
waves  in  intermediate  depth,  which  extends  the  range  of  applicability  of  the 
mild-slope  approximation  by  providing  for  relatively  rapid  undulations  in 
depth.  Deviations  from  the  mean,  slowly  varying  depth  must  be  small  in 
relative  amplitude  but  may  be  of  any  arbitrary  form.  The  present  results  thus 
extend  the  previous  analytic  results  for  sinusoidal  topography,  and  make  it 
possible  to  directly  handle  physically  realistic,  one-  or  two  dimensional  bed- 
forms. 

Although  the  results  of  Section  3  show  that  the  small-amplitude  theory  is 
able  to  predict  physically  realistic  results  for  wave  reflection  over  bed 
undulations  with  amplitudes  as  large  as  32  percent  of  the  mean  depth,  it  is 
expected  that  some  limitation  to  the  present  theory  would  occur  with 
increasingly  higher  bed  forms.  The  limitations  of  the  small  amplitude  theory 
have  been  investigated  by  Dalrymple  and  Kirby  (1984),  using  a  boundary 
integral  approach  for  bottom  undulations  of  arbitrary  height.  However,  the 
result3  of  Section  5  indicate  that  fairly  large  bottom  variations  may  be 
treated  with  a  reasonable  degree  of  accuracy. 


profiles  of  trapped  waves  for  cosine  topography 


The  symmetric  and  anti-symmetrlc  cases  correspond  to  even  and  odd  mode  numbers 
n,  respectively.  As  in  the  previous  example,  the  wavenumber  Jl  along  the  ridge 
must  be  >  kj  away  from  the  ridge,  while  the  presence  of  regions  of  constant 
depth  to  either  side  of  the  ridge  implies  the  existence  of  cutoff  conditions 
for  the  highest  mode  number. 

Numerical  results  were  obtained  for  the  particular  case  of  D/hj  =  0.3. 

The  curves  representing  the  dispersion  relation  are  given  in  Figure  II. 9. 

These  results  were  obtained  by  varying  L  and  using  a  fixed  depth  hj  =  1,  for 
which  kj  =  1.2051.  As  expected,  the  number  of  trapped  modes  increases  with 
increasing  L  (increasing  £*,  above).  The  present  computations  are  for  a 
fairly  large  amplitude  bottom  deviation;  decreasing  the  amplitude  D  would  have 
the  effect  of  decreasing  the  possible  number  of  trapped  modes  at  a  fixed  value 
of  L,  due  to  the  concurrent  reduction  in  the  £*  parameter. 

Plots  of  the  surface  profiles  for  symmetric  modes  and  anti-symmetric 
modes  are  given  in  Figures  II.  10a  and  II. 10b,  respectively,  for  the  case  of 
k^L  =  20.  In  contrast  to  the  results  of  Mclver  and  Evans  (1984),  where  the 
depth  at  large  x  goes  to  infinity,  it  is  seen  that  the  maximum  amplitude  along 
a  given  mode  profile  does  not  necessarily  occur  at  the  antinode  in  shallowest 
water,  but  instead  typically  occurs  at  an  antinode  nearer  to  the  location  of 
the  turning  point  for  the  profile  in  question.  This  effect  was  also  noted  by 
Dolan  (1983),  who  calculated  edge  wave  modes  over  real  profiles  taken  from  the 
Chesapeake  Bay.  The  profiles  studied  by  Dolan  typically  had  fairly  steep 
foreshores  and  then  were  extremely  flat  out  to  the  region  of  the  caustic 
(turning  point)  for  the  edge  wave  mode.  In  both  that  case  and  the  pr  nt 
case,  the  effect  of  shoaling  between  the  line  of  symmetry  (or  shoreline)  and 
the  region  of  the  caustic  Is  fairly  mild,  and  the  maximum  amplification  occurs 
near  the  caustic  due  to  the  refractive  singularity.  These  effects  are  in 


equivalently,  0  <  D/hj  <  0.4.  This  result  Indicates  that  the  present  small-* 
amplitude  bottom  variation  theory  Is  reasonably  accurate  In  comparison  to  the 
mi ld-slope  theory  for  bottom  variations  on  the  order  of  0-40%  of  the  mean 
depth,  even  when  the  deviation  6  Is  nearly  constant  and  extends  over  a  large 
physical  distance.  It  Is  therefore  understandable  that  agreement  between  data 
and  theory  was  good  for  the  D/h  =  0.32  case  In  Section  3. 

As  a  final  example,  we  consider  the  waves  trapped  over  a  ridge  with 
continuous  profile 


Dcos(hx/2L)  |x|<L 


S(x)  = 


0 


M>L 


(5.15) 


and  with  h  =  constant.  The  dispersion  relation  and  surface  profiles  are 
obtained  numerically  using  (5.2)  and  the  formulation  of  Kirby  et  al,  without 
any  coordinate  stretching.  Due  to  symmetry,  the  solution  is  obtained  only 
over  the  region  0  <  x  <  xM,  with  the  boundary  condition 

♦(xM)  =  0  (5.16) 

Imposed  as  an  approximation  of  the  requirement  that  M  0  as  x  +  ■".  The 
resulting  trapped  waves  may  be  symmetric  or  antisymmetric  about  x  =  0;  the 
approximate  boundary  conditions  are: 


'J'(O)  =  0  ;  anti-symmetric  case 


|>x(0)  =  0  ;  symmetric  case 


(5.17) 
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III. 2. 2  Explicit  Results  for  Partial  Standing  Waves 


We  now  wish  to  construct  parabolic  approximations  for  the  amplitude  of 
two  wave  components  of  the  same  frequency  u)  and  with  an  angle  of  180°  between 
the  assumed  direction  of  the  wavenumber  vectors.  We  choose 


=  ris. 

2u> 


A(ux,uy,ut) 


+* 


zM. 

2co 


B(px,uy,ut) 


E  + 


* 


(2.16) 


where 


_  ei(kx-ut)  .  E  =  el(-kx-wt) 


(2.17) 


and  where  *  denotes  the  complex  conjugate  of  the  preceedlng  term.  The 
propagation  direction  of  <|>+  is  oriented  with  the  +x  direction.  Using  this 
splitting  and  the  equations  of  the  preceedlng  subsection,  we  obtain 


n i  =  y  E+  +  *  +  ~  E_  +  *  +  0(p ) 


(2.18) 


-3iw  2  2  -3iu  22  .  , 

>2  =TTA  E+  +  *  ~Te~B  E-  +  *  +  0(p) 


(2.19) 


I'  ' 

k  (cosh2kh  -2)  f  ,2^2  ,  „2„2  ,  2  r 

°2  -  4  — Lnh  2kh  ~ '  tA  E+  +  B  E-  +  *•  '~*2. 


k  coshkh(cosh2kh+2'  , . 2„2  L  „2„2  ^ 

- — - -  |A  E+  +  *  B  E_  +  *)  +  0(u  ) 


8 


slnh  kh 


(2.20) 


*2 -=*■♦! <“**.**-**>  <2-2l> 


The  third  term  in  (2.21)  corresponds  to  a  rapid  standing  variation  at  cos(2kx) 
corresponding  to  the  envelope  of  the  partial  standing  wave.  Farther, 


y  ?  = 


-gh(2cosh2kh-l) 
2  sinh  2kh 


{ABE+E_  +  *} 


(2.22) 


which  represents  a  spatially  slowly  varying  oscillation  at  cos(2at).  We 

remark  that  Y2  may  not  be  eliminated  by  choice  of  b2,  li  contrast  to  the 

progressive  wave  problem.  Elimination  of  (n2  +  b2)  between  (2.13)  and  (2.14) 

2 

yields  the  forced  wave  equation  for  the  0(e  )  mean  motion: 


n  t  -v(*wh*2>  -  Sr  (l<i4i"  -  +  l»i‘>t 


gkcosh  2kh  , 

2sTn'h  ~2kh~  <AB*V-  +  *>, 


(2.23) 


The  third  terra  in  the  r.h.s.  of  (2.23)  indicates  that  spatially  fast  (~2kx) 

adjustments  in  b2  and  <)>  2  will  occur  in  response  to  slow  temporal  changes  in 
4*  — 

the  amplitude  of  4>  ^  or  <j> 

Substitution  of  (2.16)  -  (2.21)  in  the  expression  for  N.L.T.  yields  the 
result 


N.L.T.  -  -  2j  {kj> ' - —  2—  ♦  '  } 

x  2wcosh  kh  t 

.  2 

+  { k$  '  +  - 2 —  '  }  <t> . 

x  2uicosh  kh  t 

-  wVDjdAl2^  +  |  B  |  24>  ^ ) 

-  Jj2k2D2(|Bj2.^  +  1  A  1  24> " ) 


3 

+  terms  proport  Iona  1  to  E 


(2.24) 


In  the  following,  a  Fourier  decomposition  will  eliminate  terms  aE  ,  which 


would  have  been  taken  care  of  by  Inclusion  of  4> ^ .  The  coefficients  and 
are  given  by 


D,  = 


cosh  4kh  +  8-2tanh  kh 
4 

8sinh  kh 


Do  = 


-(16  sinh^kh  +  4  tanh^kh) 
8sinh^kh 


(2.25) 


(2.26) 


and  have  the  deepwater  asymptotes  (kh+")  of  1  and  -2,  respectively,  In 

agreement  with  the  results  of  Benney  (1962)  and  Roskes  (1976).  Noting  that 

the  l.h.s.  of  (2.15)  Is  a  linear  operator  on  <^,  we  then  arrive  at  decoupled 
+  — 

equations  for  <p  ^  and  4*  ^  given  by 


4;  -  v.(cc  v^>  +  (w  -k  cc  )  ♦;  +  2w  ik^  -  'aJ  ~3h  ziTh  *2  »  ♦ 

tt  x  t 


,+ 

1 


+  w2k2(D1|A|2+  D2  | Q | 2 )  =  0 


(2.27a) 


■7  -  7-(CC  7*7)  +  (w2-k2cc  )  ♦;  -  2w  {k*I  +  - ♦  '  \  4, 

®  ®  x  2wcosh  kh  t 


tt 


+  w2k2(D1 1  B  |  2+  D2  |  A |  2 )  =  0 


(2.27b) 


The  change  In  sign  of  the  terms  M>  ’  is  related  to  the  different  directions  of 
.  _  2x 

propagation  of  '{>  j  and  4>  j •  Since  D2  is  always  <0,  while  Dj  is  >0,  it  is 
apparent  that  the  presence  of  a  reflected  wave  component  <f>i  weakens  the  effect 
of  nonlinear  self-interaction  in  the  incident  wave  component.  Parabolic 
approximations  governing  A  and  B  follow  by  neglecting  time  dependence  for 


mt  1 
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purely  harmonic  motion,  and  further  by  assuming  that  0( |AXX| )<<0(k|Ax| )  as  in 
Kirby  and  Dalrymple  (1983).  The  decoupled  parabolic  equations  are  then  given 
by  (following  the  method  of  Kirby  and  Dalrymple  (1983)) 

2ik.CC  A  +  2k(k-kA)CC  A  +  i(kCC  )  A  + 
g  x  07  g  g7*^ 

+  (CCgAy)y  ‘  ^OlM2  +  °2  I B I  2  +  ^  ^Px)  A=0  (2. 28 

and 

2ikCC„B  -  2k(k-kn)CC  B  +  i(kCC  )  B  - 

o  *  u  %  ft  x 

-(CCgBy)y  +  u.VtDjBj2  +  D2|A|  2  -  ^  (4»*)x)  B=0  (2.29 

where  kQ  is  a  constant  wavenumber  defining  the  reference  phase  function  kgX  - 

oot.  These  equations  are  sufficient  for  modelling  the  propagation  of  the 

+  — 

individual  components  <j>  ^  and  <j>  ^  including  their  nonlinear  interaction; 
however,  we  have  not  yet  achieved  a  coupling  capable  of  predicting  the 
appearance  of  a  gradually  reflected  component  <j>  ^ ♦  This  coupling  is  achieved 
in  the  following  section. 

III. 2. 3  Kffect  of  Reflected  Wave  and  Mass  Transport  on  the  Nonlinear 
Dispersion  of  the  Incident  Wave 

The  effect  of  reflected  wave  components  and  induced  flow  on  the  third 
order  dispersion  of  the  incident  wave  may  be  examined  explicitly  in  the  case 
of  one  directional  propagation  over  topography  varying  only  in  the  direction 
of  propagation.  In  this  case,  we  may  integrate  (2.23)  after  neglecting  time 


dependence  to  obtain 


(2.30) 


♦  J  "  Sh  (  lA|2  "  lB|2)  +  Constant 
x 

where  the  Integration  constant  represents  an  externally  Imposed  flow  which  is 
set  to  zero.  Assuming  that  j B  j  <  |  A | ;  i.e.f  that  <j>  ^  represents  a  small  wave 
arising  due  to  reflection,  and  taking  R  =  j  B  J / J  A  j  to  represent  a  reflection 
coefficient,  the  governing  equation  for  the  incident  wave  may  be  written  as 

2tkCCgAx  +  i(kCCg)xA  -  to  2k2D*  J  A |  2 A  =  0  (2.31) 

where 

D*  -  D,  +  D2r2  ~  O  -  R2)  (2.32) 

to  h 

The  last  terra  in  D*  represents  the  effect  of  wave-induced  return  flow 
balancing  the  shoreward  mass  flux  of  the  waves.  Since  D£  <  0  for  all  kh,  the 
generation  of  a  reflected  component  reduces  the  dispersive  effect  of 
nonlinearity.  Conversely,  incrasing  reflection  tends  to  reduce  the  effect  of 
the  wave-induced  return  current,  as  the  mass  flux  of  the  reflected  wave  tends 
to  balance  out  the  flux  of  the  incident  wave. 

’k 

Figure  III.  I  gives  plots  of  D  (kh,  R)  for  a  range  0  ^R  <_  1,  with  and 

without  current  effects  included.  It  is  noted  that,  for  all  reasonably  small 

•k 

values  of  kh,  the  mean  current  has  a  noticable  effect  on  D  ,  with  the  effect 
becoming  pronounced  for  kh  <  ~  4.  It  is  also  apparent  that  the  reflected  wave 
can  act  to  entire  cancel  the  amplitude  dispersion  of  the  incident  wave,  with 
the  critical  value  being  given  by 


4^P 


T 


Neglecting  current  effects,  this  becomes 


-1/2 

which  Ls  defined  for  all  kh  and  +2  as  kh  When  current  Is  Included, 

Rc  Is  undefined  for  kh  0.649;  the  deepwater  asymptotic  value  ls 
unaffected.  A  plot  of  Rc  with  and  without  currents  Is  given  in  Figure 
III. 2.  Since  it  is  apparent  that  the  effect  of  the  wave  induced  flow  can 
significantly  alter  the  coefficient  of  the  nonlinear  term  for  values  of  kh 
small  enough  to  give  significant  depth  change  effects  over  a  short  spatial 
scale,  we  investigate  the  effect  of  retaining  or  dropping  the  current  in  a 
one-dimenslonal  example  in  Section  4,  eliminating  <j>2x  as  was  done  here.  The 
general  problem  of  coupling  between  A,  B  and  $2^  in  two  dimensions  is  left  to 
a  subsequent  study. 

III. 2.4  Extension  to  the  Case  of  Rapid  Bed  Undulations 

The  theory  derived  to  this  point  is  limited  in  *L  plication  to  waves 

propagating  over  mild  bed  slopes,  with  slope  parameter  y  presumed  to  be  of 
2 

0 (e  )  in  order  to  separate  nonlinear  and  bottom  slope  effects.  Reflections 
from  such  mild  slopes  are  likely  to  be  extremely  weak  in  intermediate  water 
depth.  Therefore,  it  is  desirable  to  incorporate  the  rapid  undulation,  small 
amplitude  bed  features  in  the  manner  of  Chapter  II.  We  thus  add  a  bed 
amplitude  parameter  k5  to  the  list  of  0(y )  parameters.  If  we  retain  the 
scaling  assumptions  used  to  this  point  in  Chapter  III  terms  like 


<5A,  <5XA,  etc. 


-  -•-  w/v 


>  v.  '- 


V  1 


.'■.'V-'y.’-V-’V: 
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Iv>  ..  . 


will  still  be  of  0(e  ;;  the  addition  of  this  feature  will  not  affect  the 

amplitude  dispersion  and  interaction  with  the  wave-induced  mean  flow  at  the 

3 

highest  order  of  magnitude  (s  )  retained  here.  We  therefore  consider  the 
extension  of  the  linear  wave  Lagrangian  formulation  to  incorporate  the 
features  of  Chapter  II.  The  nonlinear  terms  are  unaffected  and  may  be 
retained  intact  from  Section  III. 2. 2. 

Defining  the  total  water  depth  as  in  Chapter  II,  we  have 


h'  =  h  -  5 


(2.35) 


where  h  varies  slowly  with  length  scale  p ,  and  where  6  may  vary  rapidly  but 
has  amplitude  scale  u.  Following  (2.1),  the  Lagrangian  for  the  inviscid 
motion  may  be  written  as 


L  =  /  p  {*t  +  y  (Vh$)2  +  y  (*z)2  +  gz}  dz 

— h  9 


(2.36) 


We  now  expand  the  Lagrangian  about  z  =  -  h  and  retain  the  0(e  )  contributing 
to  the  linear  motion. 


Cl 


2  0  <V> 

L  =^-+r'*t  _  +  /  “1 — dz 

z=0  -h 


(V) 


o  (4>  )  (C)2 

+  /  -f—  dz  -  — 6 

-h  -h 


(2.37) 


We  now  introduce 


<P  =  f  ^Px.z)  $  t(x.i) 


(2. 38) 


‘ -V-  ’  •  -  ,  •  .■<  .  -  .  -  ••  A  *.  A'.  a'.VVV'.' 


where  f^  is  given  in  (2.5).  Making  the  substitutions  and  integrating  gives 


CC  O?j02  w2-  k2CC  r 2 


g  '  2 


6  (Vh*) 

2  2 
cosh  kh 


(2.39) 


Taking  variations  and  eliminating  n  leads  to  the  modified  wave  equation 
(II. 2. 11), 


>  -  v  .(CC  V J)  +  (u)2-  k2CC  )  $  +  - ^ - v  .(57  *)  =  0(u2)  (2.40) 

tt  h  g  h  g  .2.  h  h 

cosh  kh 


The  nonlinear  model  follows,  to  0(e),  by  adding  the  N.L.T.  terms  (2.24)  to 
the  modified  equation  (2.40). 


III. 3  The  Coupled  Parabolic  Equations 

The  desired  set  of  coupled  equations  is  derivable  in  a  large  number  of 
ways  (see  McDaniel  (1975),  Corones  (1975),  Radder  (1979)  and  Liu  \nd  Tsay 
(1983)  for  examples)  involving  the  application  of  a  splitting  matrix  to  the 
second  order  wave  equation.  We  retain  the  added  term  for  bed  undulations 
derived  above  for  generality.  Restricting  attention  to  harmont  -  waves  and 
substituting 


(1.1) 


into  (2.  15),  we  write 


(*>xx  +  p"lpx  (<nx  +  ^  =  V  +  V' 


0.2) 


-  F  —  a  ^  •  -  %  -  •  .  .  .  «  *  «_•  .  ^  '  -  V  >  *  *.1  mJ" 


where,  to  leading  order 


P  -  CC2 


'<$ 


(3.3a) 


=  (CC  )~l  [  1  +  4  (— )  6]  +  0(u  )  , 


(3.3b) 


ft'  =  gk/4ojcoshzkh  , 


(3.3c) 


,2*  -  *2  it  i  *  4(|^) «) ;  ♦  jj-  <ccg;y)y  -  ij  &  <Hy)y , 


k  g 


'l  *  55t-CD|lAl2  ^  d2!bI2  - 

g 


(3.3d) 

(3.4a) 


s2  •  6r1(»1i»i2  *  d2w2  -  k  <*Px) 

g 


(3.4b) 


We  remark  that  the  pseudo-operator  Y  is  constructed  without  isolating  a 
refraction  factor  k/k~  alleviating  the  need  for  making  the  assumption  that 


1  -  eg-)  «  1 

K0 


(3.3) 


when  performing  the  binomial  expansion  of  the  pseudo-operator  (Liu  and  Tsay, 
1983).  This  assumption  may  be  violate'  drastically  for  waves  propagating 


through  regions  with  large  depth  variations. 

We  follow  the  heuristic  scheme  of  Chapter  II  for  obtaining  the  parabolic 
equations.  Let 


>  =  iY$  +  F  +  aN^ 


(3.6a) 


<t>  »  -1y  $  -  F  +  3  N  J 

x  2 


(3.6b) 


where  F,  a,  and  0  are  undetermined,  ct  and  B  are  chosen  so  as  to  eliminate 
nonlinear  terras  from  the  coupling  term  F,  leading  to  equations  of  the  form  of 
(2.28-29).  Repeated  substitution  of  (3.6)  In  (3.2)  then  leads  to 

a  =  -  8  =  -i/2k  (3.7a) 


and 


(kCC  )  ot  ...  . 

f  <♦  -  ♦*> 
g  g 


(3.7b) 


where  we  have  retained  terms  only  to  0(p,e)  and  have  used  the  fact  that 


CCg  »  N2*  K  =  lkCCg  '  N2*  ^  +  0(G  ,y) 


(3.8a) 


and 


(C(Vx  tv*  » 


(3.8b) 


The  resulting  coupled  parabolic  equations  are  then  given  by 


it  W  -  dbr  (cciD.  *  G*  -  il  -  ™  (£)  G*  *  £  (V-)  («0 


2kCC 


gVy  2kcc 


g 


k  '•C 


y'y 


~  &  6x^+  ‘  ♦  )  +  k  V+  =  0 

g 

and 


(3.9) 


’x  +  tk*  +  2kcc  ^YyJy  2kcc 


fee  <(•")  - 

'•  Cl  \T'  \J 


<kCCf^  (♦*-♦-).  ak  (I1)  «♦'  -  rtl1)  («\)„ 

g  g  g 


Of  'x 

>  (a-i  - 

g 


)  2k  N2^  =  ° 


(3.10) 


F.quatlons  of  the  form  (2.28-29)  are  recovered  by  making  the  substitutions 


~  + 
♦ 


Ae 


lk0x 


(3.  11) 


yle Iding 


4jO  » 

2ikCC  A  +  2k(k-k0)  CC  A  +  (CC  A  )  +  2 ufl  •  [  2k<S  -  id  ]A  -  -  (6A  ) 

§x  US  S  y  y  x  k  y  y 


i (kCC  )  A  -  CC  N.A  =  { i (kCC  )  -  2iuiVd  }  Be 

g'x  g  1  g'x  x; 


-2ikQx 


(3.12) 


nd 


2ikCCgBx  -  2k(k-kn)  CCaB  -  (CCCTBV)V  -  2uSl  ’  [ 2kd  + 


^0 


g  y'y 


id  ]  B  + 
x J 


ImQ.  ' 


(dB  ) 

y  y 


2ikf,x 

+  i(kCCg)x  B  +  CCgN2B  =  {i(kCCg)x  -  2io>n'dx}  Ae  U  (3.13) 

iquations  (3.12-13)  may  be  used  in  an  iterative  fashion  to  calculate  the 
avolution  of  the  amplitude  envelopes  A  and  B.  The  numerical  scheme  used  in 
subsequent  sections  is  based  on  the  Crank-Nicolson  method,  with  each  equation 
being  solved  for  the  entire  domain  using  the  scheme  of  Yue  and  Mei  (1980), 
after  which  iteration  between  the  equations  is  performed  according  to  the 
method  provided  by  Liu  and  Tsay  (1983).  Details  are  thus  omitted. 

III. 4  Effect  of  Mass  Transport  Terms  on  Nonlinear  Reflection:  Normal 
Incidence  on  1-D  Topography 

As  a  test  of  the  effect  of  nonlinearity  on  the  partial  reflection 
process,  we  first  study  the  reflection  of  normally  incident  waves  propagating 
over  a  continuous,  one-dimensional  region  of  slowly  varying  depth.  This 
reduction  of  the  problem  allows  for  the  direct  integration  of  the  forced  wave 
equation  (2.23)  as  given  by  (2.30).  Neglect  of  the  Integration  constant 
forces  the  wave  induced  return  flow  to  balance  the  net  mass  transport,  and  is 


thus  consistent  with  the  situation  of  wave-tank  experiments.  Further,  w 

remark  that  the  generation  of  free-long  waves  would  be  absent  from  the  problem 
due  to  the  lack  of  wave-like  modulations  in  the  amplitude  envelopes. 

The  choice  of  computational  examples  which  exhibit  significant  reflection 
over  a  short  spatial  scale  is  difficult  in  this  problem,  since  water  depth 
would  have  to  be  relatively  shallow  in  order  for  an  isolated  topographic 
variation  with  an  extent  of  one  to  two  wavelengths  to  have  a  significant 
effect  on  the  incident  wave.  For  example,  the  numerical  example  chosen  by  Liu 
and  Tsay  (1983)  was  run  for  a  water  depth  at  infinity  corresponding  to 
kh=0.42.  The  range  of  admissible  waveheights  yielding  an  Ursell  parameter  of 
a  permissible  size  for  the  Stokes  theory  to  be  valid  is  thus  severely 
constrained.  Likewise,  the  linear  transition  studies  by  Booij  (1983)  covered 
a  range  of  depths  corresponding  to  0.2<kh<0.6,  again  too  small  to  be  of 
particular  use. 

For  the  purposes  of  this  study  we  have  chosen  to  consider  the  case  of 
reflection  from  an  isolated  patch  of  sinusoidally  undulating  topography  as 
studied  recently  by  Davies  and  Heathershaw  (1984)  and  Mei  (1985)  as  well  as  in 
Chapter  II. 

The  topography  is  given  by 


'  4 


h’(x)  =  1 

l 


h  ;  x<0 

hj-Dsin(2itx/ i) 
h  i  ;  x>n  l 


0<x<n£. 


(4.1) 


where  l  is  the  ripple  length  and  n  the  number  of  ripples.  The  topography  is 
illustrated  in  Figure  III. 3. 
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dimensional  patch  of  ripples  of  finite  extent  in  the  x  and  y  directions. 

Ripples  with  length  Z  are  aligned  with  crests  parallel  to  the  y-axis.  The 

patch  is  symmetric  about  the  x-axis  with  dimensions  nlx2nZ  in  x  and  y,  where  n 

is  the  number  of  ripples  as  before.  The  topography  is  specified  according  to 

Section  II. 4,  equation  (4.16).  Computations  were  run  with  2k/X  =  1,  D/h^  = 

0.3  and  a  farfield  relative  depth  kh^  =  1,  giving  h^  =  it  *  with  Z  =  1.  The 

2 

Ursell  number  Ur  =  (a/hj)/(khj)  and  wave  steepness  e  =  ka  are  thus  both  given 
by  xAq  for  the  incident  wave. 

Results  were  shown  for  the  case  of  n  =  4  in  Figures  III. 7  and  III. 8. 
Figure  III. 7  gives  results  for  the  linear  case,  with  normalized  amplitude 
I a/a0 |  given  in  7a,  |b/Aq|  in  7b,  and  the  total  wave  field  in  7c.  In  Figure 
III. 8,  results  are  shown  for  the  case  of  e  =  Ur  =  0.2.  A  comparison  of  the 
figures  indicate  some  differences  in  the  transmitted  wave  field  over  and 
downwave  of  the  ripple  patch,  with  the  amplitude  downwave  of  the  last  ripple 
being  increased  in  the  nonlinear  case,  indicating  the  greater  tendency  towards 
diffraction  effects  due  to  nonlinearity.  This  result  is  consistent  with  the 
phenomenon  of  self-defocussing  as  demonstrated  by  Kirby  and  Dalrymple  (1983), 
and  further  is  in  agreement  with  Yue’s  (1980)  results  showing  that  diffraction 
of  waves  into  a  shadowed  region  proceeds  more  quickly  in  the  nonlinear  case. 
The  reflected  wave  amplitudes  are  quite  similar  for  both  cases,  with  a  minor 
increase  in  peak  amplitude  at  the  upwave  end  of  the  ripple  patch  being  noted 
in  the  nonlinear  case. 

Comparisons  of  normalized  amplitude  for  the  total  wave  field  along  y  =  0 
and  along  x  =  31  are  given  In  Figures  III. 9  a  and  b,  respectively,  for  the 
linear  and  nonlinear  cases.  Differences  are  largely  confined  to  a  reduction 
in  transmitted  wave  height  downwave  of  the  ripple  patch  in  the  nonlinear  case. 


as  discussed  above 


(4.8) 


2ga  (kU  , 

9^  *  coskx  cos  x  +  u)t } 

g 

which  reduces  to  the  usual  stationary  form  at  U+0.  It  is  clear  that  the 
stationary  component  cos(kx)  is  unaffected  by  the  presence  of  U;  the  relation 
between  the  wave  envelope  and  topography  is  thus  unaltered. 

The  effect  of  nonlinearity  is  thus  limited  to  the  lengthening  of  the 
incident  wave  (for  small  R)  with  respect  to  the  topography  at  fixed  m 
(or  k).  For  large  2k/X,  the  effect  is  a  downshift  of  the  maxima  and  minima  of 
R  which  increases  with  increasing  e.  It  is  remarked  that,  for  large  R,  the 
incident  waves  may  be  shortened  by  nonlinearity;  the  shift  in  the  pattern  of  R 
would  then  be  expected  to  be  towards  higher  values  of  2k/X. 

In  the  present  example,  the  large  reflection  coefficient  at  the  peak  2k/X 
causes  a  significant  weakening  of  the  nonlinear  dispersion  in  the  incident 
wave  on  the  upwave  side  of  the  ripple  patch.  Referring  to  Figure  III.l,  we 
see  that  a  reflection  coefficient  of  0.6-0. 7  would  lead  to  an  effective 
nonlinear  parameter  D*  with  values  in  the  neighborhood  of  0.  Nonlinearity 
thus  has  little  effect  on  the  reflection  process  when  reflection  is  strong; 
this  is  born  out  by  the  result  of  little  or  no  shift  in  peak  value  or  location 
of  R  for  this  example.  This  conclusion  in  partially  supported  by  the 
experimental  results  of  Davies  and  Heathershaw,  who  saw  little  or  no  effect  of 
varying  wave  steepness  on  peak  reflection,  up  to  the  point  of  breaking  in  the 
incident  wave. 

III. 5  Two-dimensional  Topography 

In  order  to  take  advantage  of  the  relatively  strong  reflections  caused  by 
the  undulating  topography  studied  in  the  previous  chapter,  we  construct  a  two- 
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The  curve  for  nonlinear  reflection  shown  In  Figure  III. 5  is  for  the  case 


of  no  wave-induced  current.  It  was  found  that  the  presence  or  absence  of  the 
wave-induced  flow  had  no  significant  effect  on  the  reflection  process;  the 
corresponding  curve  for  reflection  including  mean  flows  is  thus  not 
included.  This  result  may  be  partially  explained  by  considering  a  simplified 
set  of  equations  taken  from  (2.28  -  2.29); 
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A  =  0 


B  =0 


(4.3) 
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where  depth  is  constant  and  amplitude  dispersion  is  neglected.  The  quantity  U 
may  represent  the  wave-induced  return  flow  or  a  small  flow  (<0(p))  imposed  by 
boundary  conditions.  Oscillatory  solutions  of  constant  amplitude  are  given  by 


(A(x),  B(x)}  =  {a, b}e 


ik'x 
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where 
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For  simplicity,  b  may  be  set  equal  to  a;  the  resulting  standing  wave  may  be 
written  as 


v'.--.  .V-V-V-W  '  >  v  V  V  ■.*  *  «'.'j  ■ 


■ 


73 


In  Chapter  II,  a  numerical  solution  to  the  linearized  elliptic  problem 
was  developed  in  order  to  compare  to  the  analytic  linear  solution  of  Davies 
and  Heathershaw  (1984).  The  elliptic  solution  was  shown  to  provide  an 
accurate  reproduction  of  the  analytic  results;  comparisons  with  data  were 
given  above  in  Figures  II. 2  and  II. 3. 

The  theoretical  results  of  Davies  and  Heathershaw  were  used  as  a  check  of 
the  iterative  scheme  (3.12-13)  in  its  linearized  form.  Results  for  the  cases 
illustrated  in  Figures  II. 3  and  II. 4  were  recomputed  using  the  coupled 
parabolic  equations,  using  a  grid  spacing  Ax=£./20.  The  reflection 
coefficients  for  both  cases  are  shown  in  Figure  III. 4.  The  peak  reflection 
coefficients  R  obtained  using  the  parabolic  model  are  slightly  less  than  the 
values  predicted  by  the  corrected  theory  of  Davies  and  Heathershaw.  For 
2k/A>2.0  the  number  of  points  per  incident  wavelength  drops  below  20,  and  some 
reduction  in  the  reflection  coefficient  was  noted.  A  further  reduction  in 
grid  spacing  to  Ax=£/10  caused  a  reduction  of  the  reflected  wave  amplitude  at 
2k/\=l  of  6.6%  for  the  case  with  n=10.  Therefore,  the  value  of  Ax=2./20  was 
used  for  all  subsequent  runs. 

We  now  use  the  nonlinear  form  of  the  parabolic  equation  to  study  the 
reflection  process.  Tests  were  performed  using  the  geometry  of  Figure  III. 4a 
(D/hj  =  0.16);  results  are  presented  in  Figure  III. 5  for  incident  wave 

2 

steepness  e  =  kAg  =  0.2.  A  plot  of  Ursell  number  Ur  =  (Ag/h)/(kh)  for  this 
case  is  given  in  Figure  III. 6.  The  plot  indicates  that,  for  the  chosen  value 
of  e,  results  of  the  Stokes  wave  model  are  only  roughly  valid  for  the  region 
2k/A<l.  For  small  values  of  2k/A  (large  wavelength)  it  is  anticipated  that 
the  incident  wave  phase  speed  is  overestimated.  This  would  have  the  effect  of 
over  emphasizing  the  differences  between  the  linear  and  nonlinear  reflection 


The  effect  of  diffraction  in  the  present  example  may  be  estimated 
simply.  If  we  take  the  peak  reflection  along  the  centerline  to  be 
approximately  R  =  0.6  from  Figures  7  or  8  b),  we  conclude  that  the  transmitted 
amplitude  downwave  of  the  patch  should  be 

2  1/2 

T  =  Cl-R  )  =  0.8 

This  value  is  in  rough  agreement  with  the  linear  result  but  underestimates  the 
nonlinear  value  of  T  =  0.9,  again  indicating  the  more  rapid  effects  of 
diffraction  in  the  nonlinear  case. 

III. 6  Discussion 

In  this  study  we  have  utilized  a  variational  principal  to  develop  a  wave 
equation  governing  the  propagation  of  Stokes  waves  in  a  varying  domain,  after 
which  use  is  made  of  a  splitting  method  to  provide  coupled  equations  for 
forward  and  back-scattered  components  of  an  initially  plane  incident  wave 
propagating  over  uneven  topography. 

The  restriction  to  Stokes  waves  and  the  resulting  constraints  on  water 
depth  relative  to  the  incident  wavelength  made  it  difficult  to  develop 
computational  examples  which  describe  a  significant  reflection  process  arising 
over  a  short  spatial  scale.  Under  the  mild-slope  conditions,  it  is  likely 
that  the  gradual  reflection  process  would  be  apparent  over  only  fairly  long 
spatial  scales.  For  the  case  of  shallower  water,  the  Stokes  wave  formulation 
is  no  longer  valid,  and  recourse  must  be  made  to  appropriate  equations  such  as 
the  Boussinesq  equations.  The  parabolic  approximation  for  a  spectrum  of 
steady  (in  time)  forward-scattered  waves  in  the  shallow  water  regime  has  been 
provided  by  Liu,  Yoon  and  Kirby  (1985);  the  development  of  a  model  for  partial 


reflection  in  this  case  will  be  the  subject  of  a  further  investigation.  Also 
of  special  interest  is  the  case  where  the  incident  wave  amplitude  is  modulated 
in  space  and  time.  The  treatment  of  "groupy"  waves  is  not  approachable  using 
the  reduced  wave  equation  of  Section  III. 3;  however,  the  general  time 
dependent  model  (2.15)  may  form  the  basis  of  such  an  approach,  after  further 
accounting  for  terms  arising  due  to  possible  fast  modulations  of  0(e)  in  the 


amplitude  functions. 


The  integrals  I  are  defined  as  the  integral  over  total  depth  of  an 
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Appendix  III.B  Components  of  the  primitive  Lagranglan  L 


The  Lagranglan  L  is  expanded  as  a  series  in  powers  of  the  wave  steepness 
parameter  e  without  regard  to  the  relative  size  of  the  modulation  parameter 
p.  After  expanding  n  and  the  integrals  I  in  (2.6),  the  individual  components 
of  L  in  (2.8)  are  given  by  (after  dividing  out  the  constant  density  p) 
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Appendix  III.C.  General  Form  for  0(e)  term  in  wave  equation 


The  term  {N.L.T.}  In  (2.15)  Is  given  here  in  terms  of  (h2  +  *>2' 
and  4>2 •  The  first  order  surface  rij  has  been  eliminated  through  use  of  (2.9) 
Further,  we  have  made  use  of  the  fact  that 
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Vh^l’^2^  =  _4k24>2)  +  0(u) 


for  both  progressive  and  standing  waves.  {N.L.T.}  is  then  given  by 
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where  we  have  substituted  for  all  I  values  from  Appendix  III. A. 
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